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We consider random generalizations of a quantum model of infinite range introduced by Emch 
and Radin. The generalizations allow a neat extension from the class h of absolutely summable 
lattice potentials to the optimal class 1% of square summable potentials first considered by Khanin 
and Sinai and generalised by van Enter and van Hemmen. The approach to equilibrium in the 
case of a Gaussian distribution is proved to be faster than for a Bernoulli distribution for both 
short-range and long-range lattice potentials. While exponential decay to equilibrium is excluded in 
the nonrandom l\ case, it is proved to occur for both short and long range potentials for Gaussian 
distributions, and for potentials of class k in the Bernoulli case. Open problems are discussed. 
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I. 



INTRODUCTION AND SUMMARY 



Equilibrium properties of general random quantum systems have been the subject of recent investigations [1]. 
Mean-field models, e.g., quantum analogues of the Sherrington-Kirkpatrick (SK) spin-glass have also been studied, 
and the thermodynamic limit proved for the SK model in a transverse external field [2]. In this paper we consider 
the dynamics of a special class of random quantum models of infinite range, which may be viewed as a version with 
random couplings of the model introduced by Emch [3] and generalized by Radin [4]. 

The most compelling reason for studying random quantum systems - in particular, quantum spin glasses - is 
that only for a quantum system does a physically satisfactory definition of the dynamics of states and observables 
exist. Kinetic theory approaches to the SK spin glass model, for instance, rely on Glauber dynamics [5], which is 
an imposed dynamics: indeed, the Ising model is, of course - both mathematically and physically - to be regarded 
as an anisotropic limit of suitable quantum systems, whose dynamics is naturally defined. The occurrence of ageing 
and other non-equilibrium phenomena were also derived under the assumption of an imposed dynamics, such as the 
Langevin dynamics for the spherical and mean spherical model with competing interactions treated in ([6], [7]). 

The Emch-Radin model seems to be largely ignored in the recent literature. Indeed, in an influential paper on the 
mathematical theory of non-equilibrium statistical mechanics [24], Jaksic and Pillet remark that, to their knowledge, 
the dynamical aspects of quantum spin systems have been little studied, the only well-understood case being the 
one-dimensional XY model [25]. An additional point lending special interest to the model of ([3], [4]) is that it has had 
a qualitatively successful application to the description of real experiments, which display a non-Markovian approach 
to equilibrium ([27], [28]), as well as the fact that several results are not restricted to dimension one ([3], [4]) . From a 
purely conceptual point of view, the quantum extension of the Kac ring model in [26] should also be mentioned. 

We now comment on the organization of the paper. 

Section II reviews the nonrandom model ([3], [4]). We show that the infinite volume limit (<rf)(t)of the one-site 
transverse spin exists if the lattice potential is square-summable (class I2). For the case in which the lattice potential 
is absolutely summable (class li), for no form of the interaction the approach to equilibrium is of exponential type: 
this is a well-known result in function theory, mentioned for completeness in proposition 2.1, because reference is 
made to it later. 

Section III introduces the random models, which furnish a natural interpretation of the result of section II for square 
integrable lattice potentials because they are, in contrast to the nonrandom case, thermodynamically stable. This 
was first observed for classical models by Khanin and Sinai [8] and extended, with an independent proof, to quantum 
systems by van Enter and van Hcmmen [9]. For definiteness, we restrict ourselves throughout the paper, to two 
opposite extreme examples of probability distributions, viz. Bernoulli and Gaussian, and show that the Bernoulli rate 
of approach to equilibrium coincides with the nonrandom case. We also prove that the mean transverse magnetization 
converges along a subsequence, almost everywhere in the probability space, to its average, a nonrandom function fit). 

In section IV we show that the approach to equilibrium is always faster than Bernoulli ,and exponential, for the 
Gaussian distribution, both for short range and long range potentials of classes l\ and I2, while for Bernoulli it is 
exponential only for potentials of class I2. 

Section V is reserved to a conclusion and open problems. 



Following Radin's review [4] of the work of Emch [3], consider the experiment of [9] (see also [10] ): a CaF 2 crystal 
is placed in a magnetic field thus determining the z-direction, and allowed to reach thermal equilibrium. A rf pulse 
is then applied which turns the net nuclear magnetization to the x direction. The magnetization in the x direction is 
then measured as a function of time. 

As in ([3], [4]), we assume an interaction of the form 



II. THE NONRANDOM MODEL 




(1) 



where j S V and k 6 V in (1), where V is a finite region, V E Z". Hy is defined on the Hilbcrt space 



(2) 
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The state representing the system after the application of the rf pulse will be assumed to be the product state 

P = p(0) = (g)^- (3a) 
jev 

where 

M-) = tr A- exp(-7^))/^(exp(- 7( Tj : ) (36) 
where tr is the trace on Cj. Other choices of the state are possible [4], but wc shall adopt (3) for definiteness. Let 

S x = l/N(V) °? (4) 

jev 

be the mean transverse magnetization, with N(V) denoting the number of sites in V. The real number 7 in (3b) may 
be chosen as in [3] such as to maximize the microscopic entropy subject to the constraint trS x p(0) equal to a constant, 
i.e., a given value of the mean transverse magnetization. Since the state (3a) is a product state, 7 is independent of 
V, and has the same value if S x is replaced by af Q , for any i 6 V. Let 

pY = U vP (0)Uv (5a) 

and 

«Mi)=pr«) (56) 

where Uy(t) = cxp(itH v ). We may write, by (5), 

K)v(t) = p(0)(U v (-t)af o U v (t)) (6) 

It is natural to define 

(af )(t) = lim (af ) v (t) (7) 

V — >oo 

provided the limit on the r.h.s. of (7) exists, as the expectation value of the local transverse spin in the (time- 
dependent) noncquilibrium state p^ of the infinite system. A weak-star limit [11] of the sequence of states Py(-) = 
tr(-U v p(Q)Uy t ) exists, by compactness, on the usual quasilocal algebra A of observables. The expectation value of 
af in the equilibrium state associated to (1) is zero by the symmetry of rotation by ir around the z axis. We may 
now pose the question whether the limit 

■^^J^ oc p^ N ^T,<) (8) 

where A denotes a finite subset of Z", which is interpreted as the mean transverse magnetization, exists. The property 
of approach to equilibrium is expressed by 

lim f(t) = (9) 

t— >oo 

Of particular interest is the rate of approach to equilibrium. In the present nonrandom case, the limit at the r.h.s. 
of (8) equals f(t) = p^icrf ), for any io, by translation invariance of p^. This is not so for random systems, in which 
case additional arguments are necessary to show the convergence of the r.h.s. of (8) for almost all configurations of 
couplings (section III). 

At this point we simplify the analysis as in [3] by restricting ourselves to the one-dimensional case. Generalization 
of some results to any dimension v is straightforward, as remarked in [3] , because they rely on the general Proposition 
1 of rcf . [4] . The Hamiltonian becomes 

n n 

H n = H v={ _ nM = 1/2 <\3 - fc IK CT fc + B E °] ( 10 ) 

j,k=—n j — — n 

where the function e(-) is assumed to satisfy: 

e(0) = (11a) 
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e(n) is monotonically decreasing in n (H^) 

and either: 

oo 

|e(n)| < oo i.e. e G l\ (12a) 

71=1 

or 

oo 

^(|e(n)|) 2 < oo i.e. e € Z 2 (126) 

n=l 

The terms containing the interaction e(j) in the operator U~ t af Q U^ l (n = V) which appear on the r.h.s. of (6) are 
of the form (see [3], pg.1200): 

[cos(e(j» - iaty^ sin(e(j)i)] x [< cos(2e(j)i) + « o+j sin(2e(j)i)] x [cos(e(j)*) + sin ( e 0»] ( 14 ) 

By (3), and emplying a product basis of eigenstates 

*a = J| With °"f *«i = a i*«i ( 15 ) 

with ojj = ±1, in which p(0) is diagonal, we see that the terms in (14) which contribute to the right hand side of (6) 
are (using a z a x — ia v et cicl.): 

((cos(e0») 2 - (sin(e(j)i)) 2 )p(0)(<)cos(2e(j)i) = = (cos(2e(j)i)) 2 «k (<) (16) 

The above formula is true, however, only if iq — j > —n, i$ + j < n in (14), due to the restriction in (10). If, e.g., 
to — j < —n, with iq + j < n, the first and third terms in (14) will be replaced by the identity, and the resulting term, 
instead of (16), will be 

cos(2e(. 7 )t)^ «) (17) 

For this reason it is important to perform the limit on the r.h.s. of (7) before the one on the r.h.s. of (8): the 
asymmetry which yields a few terms of the form (17) in the case of a finite system for io ^ disappears upon taking 
the limit n — > oo. For this reason it is also important to take the infinite volume limit in (7) in the sense of Fisher, 
which means roughly that the domains V have to increase at about the same rate in all directions (see [8] for comments 
and references): in one dimension, this condition is satisfied by the symmetric choice in (10). 

With the above remarks, we are allowed to restrict ourselves to i = in (7) and obtain, with (10), (3): 

n 

(og)n(t) = M°o) n(cos(2e(j)t)) 2 x cos(2Bt) (18) 

We now have: 

Proposition2.1 Let 

MO^o (19) 

Then, if e e l 2 , i.e., (12b) holds, 

3 lim <og) n (t) ee (ag)(t) (20) 

n^oo 

If e e ?i, i.e., if (12a) holds, for no constants < a < oo,0 < C < oo, the inequality 

K^o>(*)l < Cexp(-|a|t) if t > (21) 

is satisfied. 

Proof: By (12b) and (lib), the product on the r.h.s. of (18) converges to the infinite product 

oo 

p=H(cos(2e(j)t) 2 ) (22a) 

3=1 
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Above, p exists by the following simple bound: by (lib), for any t there exists M < oo such that 2|e(j)|i < lVj > M 
and thus 

oo 

P < d n (I cos(2i(e(j)) 2 )|) < d H (1 - ct 2 (e(j)) 2 ) < dcxp(-ct 2 ]T (e(. ? )) 2 ) (23) 

j=M 3>M 3>M 

Above, c and d denote positive constants which may vary at each step and we used the elementary inequalities: 
< cos(x) < 1 — cx 2 , 1 — x < exp(— x), for some c>0if0<x<l. We have thus proved (20), with 



(o$){t) = M°5)P (226) 

By ([4],pg.2954, and Proposition 1): 

K)(t) = (n p (og)$ p ,ex P (*tn p (2ff )* P ) (24) 

where 

oo 

Ho = *5 E <\ k ^k ( 25 ) 

k— — oo 

Above, i?o is the norm limit, as n — > oo (in the quasi-local algebra .4 [11]) of the sequence of elements 
<7 nEL-ri £ (l^l) (7 t' which is the generator of the group of time-translation automorphisms of A, and p = the 
state on A, p — ®jez0j corresponding to (3a), with H p the GNS representation associated to p with cyclic vector $ p 
[11]. By (25), Hq is bounded below (in fact it is bounded) and (21) follows from a standard result in function theory 
(see,e.g., [10], 5.8). 

Remark 2.1 

The proof of non-exponential decay relies on the form (25), which is meaningless when only (12b) (but not (12a)) 
is valid, and, indeed, we see from (23), by choosing M > t 1 / , that, for potentials of class Z 2 of the form (48a), with 
1/2 < a < 1, exponential decay (in the sense of inequality (21)) is possible when only (12b) is required. The latter 
suffices to prove (20), as we saw, but it is clear from (2) or (10) that it precludes thermodynamic stability for the 
present model. We shall see in the next section that this is remedied by the introduction of randomness. 

III. THE RANDOM MODEL 

We now introduce the Hamiltonian of a disordered system corresponding to (10) (with B = for simplicity): 

n 

H n = l/2 J(j,k)e(\j-k\)a]a z k (28) 

j,k=—n 

where J(j, k) are, again for simplicity, independent, identically distributed random variables (i.i.d. r.v.). We shall use 
Av(-) to denote averaging with respect to the random configuration J. The J are assumed to satisfy [8]: 

Av(J(j, k)) = (29a) 

\Av(J n (j,k))\ < nlc n \/n = 2,3,4,- •• (296) 
Let the free energy per site /„ be defined by 

— kT 

fn(J) = log Z n { J) (30a) 

where 

Z n (J) - tr(exp(-0H n )) (306) 
is the partition function and the trace is over the Hilbert space (2). Then 
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Theorem 3.1 [9] Under assumptions (12b) and (29), the thermodynamic limit of the free energy per site 

f(J) = lim f n (J) (31a) 

n— >oo 

exists and equals its average: 

f(J) = Av(f(J)) = lim Av(f n (J)) (316) 

n^oo 

for almost all configurations J ( a.c. J). 

In [9] general v dimensional systems were considered (as remarked, we could equally have done so here), with 
the limit of finite regions taken in Fisher's sense. The reason why (12b) suffices for the existence of the ther- 
modynamic limit is that, in order to obtain a uniform lower bound for the average free energy per site, the 
cumulant expansion (see, e.g., ([13], (12.14), pg 129, for the definition of Av c , there called Ursell functions): 
Av(exp(t J(i,j)) = cxp(^^L 2 Av c (J n (i, j))t n /n\) was used [8], which, by (29a), starts with the second cumulant 
Av c (J 2 (i, jj) = Av(J 2 (i,j)), which is the variance of J(i,j). Condition (29b) was used to control the sum in the 
exponent above. 

We now investigate whether model (28), which is thcrmodynamically stable in the sense that the thermodynamic 
limit for the free energy exists under assumption (12b) for the lattice potential by Theorem 3.1, also yields a well defined 
thermodynamic limit for the dynamics, i.e., of (crf )(t). Instead of (14), a typical term containing the interaction e(j) 
in the operator U~ t af Q U^ l appearing on the r.h.s. of (6) is, now, 

(cos(t J io - jt io e(j)) - sHJto-jMJ)*)) x « cos ( 2J io,io+je(j)t) + «+i sin(2J io;io+j e(j)t)) 

x (cos(J io _ jiio e(i)i) + sm(J io _ jiio e(j)t)) 



(32) 

We ignore the differences commented after (15), which are analogous in the present case but, again, will disappear 
upon taking the limit lim^oo as discussed there. For any fixed configuration J, an estimate of the same type (23) in 
proposition 2.1 yields a formula analogous to (22b) in the limit n — > oo: 

(a?)(t)=PU<)=M<r?)Pi (33) 

where 

oo 

^ = IT C0S (2tJi,i+k£(k)) cos(2tJi- kyi e(k)) (34) 
fe=i 

and we use the same notation as in (8). We now obtain for the finite version of the r.h.s. of (8), with A = [— m, m] 
and N(A) = 2m + 1: 



Em x sr-^m 

P^m, = Pool 2m + 1 ) = = S( 2m + 1 ) (35) 

where 

5 = U°i) (36) 

is a nonzero number independent of i. We have now: 

Theorem 3.2 Let the r.v. J be Bernoulli or Gaussian. Then there exists a subsequence ™ r= i.oc such that 



lim p^(m r ) = p io ee Av(p t0 ) (37) 



r* 



a.e. J. 

Proof. Consider the (2m + 1) x (2m + 1) matrix with the r.v. pi as entries, pi and pj, for i ^ j, arc not 
independent. For, let j > i. Then, 



pi = ]^[ cos(2t J iti+k t(k)) cos(2t Ji_ fcji e(fc)) (38a) 



fe=i 
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and 

oo 

pj = cos(2i Jjj + ke(k)) cos(2tJj-k,je(k)) (386) 
fe=i 

If i + fc = j, j — k = i, and there will be one common factor cos(2tJije(k)) in pi and pj. We now consider the r.v. 
on the r.h.s. of (35): 

X m s (39) 



Since the Jij are identically distributed, pi is independent of z. Consider 



We have: 



Xm = — 2^n — (40) 



lim Av((pi - pi)(pj - p~j)) = (41) 



To prove (41), let j > i for definitcness. By (38), with k = j — i, 

\Av((pj - pj)(p k - p k ))\ < const. \Av(cos(2J hl+k e(k)t)) 2 - (Av(cos(2J hl+k e(k)t))) 2 \ (42) 
The r.h.s. of (42) is identically zero for Bernouilli r.v.. For the Gaussian, 

A « for L°lcxp(-x 2 )(cos(2x e (fc)^) 2 l + 4exp(-t 2 ( e (fc)) 2 ) 
Av((cos(2J lyl+k e(k)ty) = ?r - = (43a) 

vW 2 

while 

Av((cos(2J^ +k e(k)t)) 2 ) = exp(-2(e(fc)) 2 t 2 ) (436) 

Putting (43) into (42) we find that the r.h.s. of (42) is 0((e(fc)) 4 ) for large k, which tends to zero by (lib) and (12b), 
thus proving (41). By (40) and (41), 

lim Av{{X m f) = (44) 

m— >oo 

Thus, by Chebyshev's inequality, X m — * in probability, which implies by the Borel-Cantelli lemma that a subsequence 
converges to zero a.e. J (see, e.g., [12], Theorem 4.2.3). This proves (37). 

The convergence a.e. J is the well-known property of self-averaging, which we conjecture holds even without the 
necessity of restricting to a subsequence. As discussed in [8] and more extensively in [13] ("Thou shalt not average") 
it is essential for the reproducibility of the outcomes of experiments, in this case the measurement of the transverse 
magnetization. 

By (8), (37) and (43b), we have the results: 

OO 

f B (t) = l[(cos(2te(k))) 2 (45a) 
fe=i 

for the Bernoulli distribution and 

oo 

/ G (t)=exp(-2£ 2 5>(£:)) 2 ) (456) 
fe=i 

for Gaussian r.v.. 

Remark3.1 By (45a) and (22a, b), the rate of approach to equilibrium agrees with the nonrandom rate in the 
Bernoulli case. This leads immediately to the interpretation of (22a, b) for lattice potentials of class I2, i.e., satisfying 
(12b), in view of Theorem 3.1. See also remark 2.1. 
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IV. THE RATES OF APPROACH TO EQUILIBRIUM 

By (45b) fa is explicit in the Gaussian case: 

oo 

/ G (^cxp(-2t 2 5>(fc)) 2 ) (46) 
fe=i 

From (46) it is exponential in the sense of inequality (21) whatever the range of the interaction. We have, now: 

Proposition 4.1 The approach to equilibrium for Gaussian r.v. is faster than for the Bernoulli distribution in the 
case of potentials of class I2, i.e.,: 

lim Ml = 00 (47) 

jG(t) 

if 

e(k) = k~ a with 1/2 < a < 1 and k > 1 (48a) 
The same holds in the prototypical example of infinite, but extreme short range case 

e(k) = 2~ k - 1 with k > 1 (486) 
Proof We first consider (48a) and use the elementary inequalities: 

cos(a;) > 1 - x 2 /2 (49a) 

valid for < x < 1 and 

exp(-ar) <l-x/2 (496) 

valid for < x < 1.5936, to find a lower bound to /s(i), given by (45a). We assume that the time values are sampled 
at points differing from the zeros of the cosine function in both cases (48a) and (48b). Then, by (49), 

00 

i/ B (t)i=in( cos ( 2 ^ Q ) 2 i^ 
fe=i 

const. {2tf' a ] j (1 - l/2(2ifc-") 2 ) 2 > const. (2i) 1 / Q exp( 8t 2 k- 2a ) 

2tfe-°<l 2tfc-°<l 

/•oo 

> const. (2i) 1/Q exp(-8t 2 / dx X - 2a ) 
= const. exp(-dt 1/Q +clog(i)) 

(50) 

where d, c are positive constants. (47) follows from (46), (48a) and (50). We now turn to (48b), for which an exact 
result exists for fs{t)'- 

00 00 
f B (t) = n(cos(2t/2 fe+1 )) 2 = [](cos(t/2 fc )) 2 = (sin(t)A) 2 (51) 
fe=i fe=i 

by Vieta's formula 



^ = g cos(t/2 fc ) (52) 



1 fe=1 

and (47) is obvious from (46) and (51). 

As a last remark concluding this section, Vieta's remarkable formula (52) may be viewed from the point of view of 
characteristic functions and hence of addition of independent r.v. in probability theory: note that the l.h.s. of (51) is 
the characteristic function of the uniform distribution in [—1, 1] ([12], exercise 8, pg.165, or [14], Chap. 1). It admits, 
however, an elementary proof, discovered by S.P.Heims([28] , Appendix): using the identity sin a = 2 sin(a/2) cos(a/2) 
successively, we obtain sin£ = 2" sin(£/2") cos(£/2" cos^/2"- 1 ) • • • cos(£/2). But lim IWOO sin ^T ] = 1, from which 
(52) results. 
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V. CONCLUSION AND OPEN PROBLEMS 

We have shown that the inclusion of randomness permits a natural generalization of the Emch-Radin model to the 
class Z 2 of square-summable lattice potentials. The generalized model (as well as the original one, of course) provides 
one of the very few examples of a random quantum system allowing a rigorous analysis of the dynamics of approach 
to equilibrium, and, in particular, of the role probability distributions play there. 

In [15], together with S.R. Salinas, we studied the phase diagram of the mean-field random Ising model, and proved 
that a tricritical point exists there for the Gaussian but not in case of a Bernoulli distribution. Thus, equilibrium 
properties are, in general, also quite sensitive to the nature of the probability distribution. We suggested [15] that the 
observed differences might be due to the fact that a discrete distribution of probabilities samples just a few values of 
the couplings and therefore introduces some short-ranged elements into the problem (sec also [16]). From this point 
of view, it may be conjectured that discrete distributions have a closer connection with real materials and, indeed, 
in the nonrandom case - which is a caricature of a real old experiment [27] - a decay of type (52) provides a better 
qualitative description than (48). The issue can only be settled by performing analogous experiments on random 
systems, as briefly discussed below. In particular, recent great progress in spin-glass theory (see [17] and references 
given there) pertains to the mean-field SK theory, which is long-range. For these, gaussian distributions may be more 
adequate, since they sample many values of the couplings, and "reinforce" in this sense the long-range nature of the 
interactions. For such systems, proposition 4.1 predicts a faster approach to equilibrium. 

Concerning the relation to experiment, the experiments at the University of Chicago and ATT Bell Laboratories 
[29] have shown the importance of the quantum regime in the study of the phase diagram of certain materials. In 
the choice of [29] , a strong spin-orbit coupling between the spins of the magnetic ion component of the material and 
the underlying crystal essentially restricts the spins to orient either parallel or antiparallel to a specific crystalline 
axis: such spins are usually referred to as Ising spins. Under application of a transverse magnetic field, oriented 
perpendicular to the preferred axis, a flipping effect of the Ising spins takes place, analogous to the one observed 
in a spin- resonance experiment, and it becomes possible to identify a quantum crossover regime at sufficiently low 
temperatures and certain (critical) values of a parameter proportional to the transverse field intensity. This has 
been nicely reviewed by S. Sachdev in [30]. The qualitative similarity to the present model is obvious, but doing an 
experiment on such materials under conditions analogous to the old experiment reported in [28] seems to be a very 
interesting open problem. 

We finally come to the interesting conceptual question of exponential versus nonexponential decay in the sense of 
inequality (21). By (48), exponential decay occurs in the Gaussian case for interactions of both short and long range, 
and in the Bernoulli case for very long-range interactions of class Z 2 by (23). On the contrary, for nonrandom systems 
exponential decay is excluded by proposition 2.1, in contrast to the behavior of the recurrence time probability in the 
classical mechanics of fully chaotic systems (see, e.g., [18], for a review and references). 

For one particle in a short-range potential in quantum mechanics, algebraic decay 0(t~ u t 2 ), where v is the dimension, 
is the rule, as known from scattering theory, due to the spreading of the free wave-packets (see, e.g., [19]). In contrast, 
algebraic decay results from mixed phase space and trapping of trajectories in regular islands [18]. However, for 
quantum systems in external (non-decaying) fields, such as constant electric fields [20] or parametric oscillators 
[21], exponential decay does occur: the Hamiltonian is not bounded below, i.e., its spectrum is the whole line and 
proposition 2.1 does not go through. Such models, dubbed quantum Anosov systems in [22], are characterized by 
positive quantum Lyapunov exponents ([20], [21], [22]). In the present model, the existence of exponential behavior is 
due to the randomness, while in the cited systems the non-semiboundedness of the Hamiltonian arises from neglecting 
the back-reaction of the field. 

A final point is that, in spite of the simplicity of the present model, it is a genuine quantum many-body problem, 
i.e., of infinite number of degrees of freedom, about which almost no results on time evolution with good ergodic 
properties exist: one exception is [23]. There it is pointed out that the condition for existence of a state-dependent 
time-evolution in the Heiscnbcrg picture is thermodynamic stability in the sense of the N-particle Hamiltonian having 
a lower bound of the form — const. N, which would require, in our case, lattice potentials of class Accordingly, 
the Heisenberg time evolution is well-known to exist as an automorphism of the quasi-local algebra in the l\ case (see 
[4] and [11]), i.e., in a state- independent sense, and it is at best state-dependent in the Z 2 case: an example is the 
present construction, which is valid for the same class of states considered in ([3], [4]). 

Acknowledgement: We thank the referees for relevant remarks and corrections. 
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